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Abstract 

Assuming a second-order phase transition for the hadronization process, we attempt 
to associate intermittency patterns in high-energy hadronic colhsions to fractal struc- 
tures in configuration space and corresponding intermittency indices to the isothermal 
critical exponent at the transition temperature. In this approach, the most general 
multidimensional intermittency pattern, associated to a second-order phase transi- 
tion of the strongly interacting system, is determined, and its relevance to present 
and future experiments is discussed. 
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1. Introduction 



One of the central issues in quark-gluon plasma physics is related to the nature 
and structure of the quark-hadron phase transition and its measurable effects on 
the pattern of multihadronic states generated in a high energy-density collision of 
relativistic ions. Recent studies in lattice gauge theories suggest that the QCD phase 
transition for two flavours (/ = 2) and massless quarks (m„ = rrid = 0) is of second 
order and this result persists even for / = 3 provided that the mass of the strange 
quark is sufficiently large [rris > 100 MeV) This development brings as to the 
hypothesis of a second-order chiral phase transition in QCD close to the real world 
and suggests a further exploitation of this fundamental assumption in an attempt to 
build-up a self-consistent dynamical scheme of strong interactions at the transition 
temperature (T = Tc). In a recent work along this line F.Wilczek has argued that 
the hypothetical second-order chiral phase transition of the / = 2 QCD belongs 
to the universality class of the n = 4 Heisenberg magnet; therefore, in order to 
study the behaviour of the strongly interacting system near the critical point, one 
may employ the critical exponents of this magnetic system as given by the fixed 
point of the corresponding renormalization group 0. In particular, the power law 
at T = Tc connecting the condensate (qq) (order parameter) with the quark mass 
(ordering field), (qq) ~ m^/^, in the limit m ^ (m = m„ = rrid) introduces a 
universal critical exponent 6 in the strongly interacting system, which refiects the 
non-perturbative sector of the underlying theory. A similar power law, (p) ~ ml/^, 
with the same exponent 6, relating the density of hadrons (pions) in configuration 
space, (p), with the pion mass (chemical potential) in the chiral limit (m7r — >■ 0), 
is expected to hold in the hadronic phase of a quark-gluon production process, at 
T = Tc. The fundamental index S has a distinct geometrical interpretation giving 
the fractal dimension. Dp = of the critical system at T = Tc, embedded in a 
support space of dimension D 0. 

In this work, adopting the hypothesis of a second-order phase transition in the 
strongly interacting system, we attempt to incorporate the above geometrical prop- 
erty of the critical system into the hadronic S-matrix for multiparticle production 
and show how the fractal structure of the critical system in configuration space may 
induce multidimensional intermittency patterns in momentum space, in accordance 
with experimental evidence. For this purpose we derive a suitable integral represen- 
tation for multiparticle densities in momentum space, introducing proper light-cone 
variables in the treatment of the reduced S'-matrix elements. In this approach the 
constraint of self-similarity (monofractality) in configuration space, required by the 
hypothesis of a second-order phase transition, is easily incorporated and its effect on 
the particle densities in momentum space is extracted. As a result, the general fractal 
solution leading to a complete multidimensional intermittency pattern in momentum 
space is obtained, expressed as a set of power laws for the factorial moments of all 
orders. This pattern depends on two universal indices, one of which is the critical 
exponent 6. Therefore, a systematic comparison with experiments may give, not 
only an indication for the onset of a second-order phase transition in the quark-gluon 
system, but also a measurement of the corresponding critical exponent at T = Tc. 

In section 2 the formalism leading to a complete set of hadronic densities in 
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3 + 1 space-time is presented by considering the inclusive S'-matrix treatment of a 
hadronic collision within the framework of inside-outside cascade dynamics [Q. In 
section 3, the geometrical constraint of fractality as a requirement of a second-order 
QCD phase transition at T = Tc is imposed and the associated fractal solution in 3D- 
momentum space is found and discussed. In section 4, the resulting multidimensional 
intermittency pattern showing strong fluctuations in small domains of momentum 
space, is derived. This pattern, combined with a smooth background representing 
conventional correlations in momentum space 0, is compared, in section 5, with 
the typical behaviour of the factorial moments observed in high-energy hadronic 
collisions. Finally, in section 6 our results are summarized and discussed and a 
critical comparison with other approaches is attempted. 



2. Hadronic densities in 3 + 1 space-time 

Consider the inclusive process 

hA + hB ^ hi + h2 + ... + hp + anything 

in the hadronization sector of the collision, defined as the component of the S'- 
matrix connecting the initial hadrons (/i^, /i^) with final states consisting of newly 
hadronized particles (at T = Tc) in the process of a second-order QCD phase transi- 
tion. The densities pdqi, in momentum space are written as follows: 

Pc(gi,.-,gp) = ^^^^"""^^ ,3^^'"^ - l(gi---gp; {^a} out|p^pB in)|^ (2.1) 

'^h Oj qi-.-CL {ka} 

where qi = {qi, Ei) are the four-momenta of hadrons (pions) hi. ..hp in the final state, 
(j/j the total cross section of the collision h^ + hs in the hadronization sector, and da^l^ 
the p-hadron inclusive cross section in the same sector. The summation runs over the 
hadronic states that build-up the missing mass in the inclusive process with a total 
four-momentum k = J2a ^a- The matrix elements in eq.(2.1) connect the initial state 
\paPb', in) with final states |out) generated at T = Tc by the quark-hadron phase 
transition (the critical component of the S'-matrix). Reducing now the S- matrix 
elements in eq.(2.1) we remove the hadrons hi. ..hp from the outgoing states |out) and 
using completeness of the states we finally obtain 0, ||: 

Pc(gi, qp) = fAB ffl d'xid''x,d^z,e-"^'''^e-''^^^^^-'^^G{xi...Xp; Z2...Zp) (2.2) 

1=2 



■ {pAPB\T{(t>{Xi)...(l){Xp)}T{J{Q)<i){z2)...<t){Zp)}\pAPB) (2.3) 

where Jab = (27r)^/(2P+^W^Pcr/i) is a normalization factor corresponding to a total 
energy W and incident momentum P = \Pa\ = \Pb\ in the centre-of-mass system. 
The current J(0) in eq.(2.3) is associated with the pion field 0(x) and denotes 
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the corresponding Klein-Gordon operator. Changing integration variables in eq.(2.2): 
x\'^^ — Xi + Zi Ti — Xi — Zi (i — 2, ...,p) and putting xi = ri, we have: 

Pc(gi, qp)-M ffl d'nd'nd'xt^e-^'^^^^e-^^^^^g{n..T,- 4+^..a;W) (2.4) 

Integrating over the variables x^'^^ we finally obtain: 

pc{qi, ■■■,%) = J d'^ri...d^rpexp{-iqiri- ... -iqprp)p{ri,..,rp) 
p{n,...,rp) = ^/nciM^^^(ri...r,;xW...4+)) (2.5) 

^ i=2 

In the hadronization sector the space-time development of the collision is assumed to 
follow inside-outside cascade dynamics; it is therefore suitable to use, in the longitu- 
dinal space-time {r\\,t), light-cone variables: r = (t^ — r\\'^y^'^,^ = |ln|^:||i^|. With 
these new coordinates we have: 

q-r^Et-q-f = ]^{E ^ q\0{t - r\{) + ]^{E - q{){t ^ r{) - q^ ■ ri_ 

ex^{iqr) — exp(— • fj,) exp[imTTCOsh(^ — y)] (2.6) 

where ttit (transverse mass) and y (rapidity) are given by the well-known relations : 
= {q\ + m^)^/^, y — The Fourier integral (2.5) is now written as 

follows: 

pc{qi, gp) = / n ndndiid^fu exp[iq±i-f±i-iTimTi cosh(yj -^j)]p(ri, Vp) (2.7) 

i=i 

where integration over the impact parameter space (r^) has been explicitly intro- 
duced. The physical picture we have adopted in this approach (inside-outside cas- 
cade mechanism and a second-order phase transition in thermal equilibrium) requires 
strong constraints on the integrand p(ri, ...,rp) in eq.(2.7). In fact since the develop- 
ment of the strongly interacting system in the longitudinal space-time is assumed to 
take place along the hyperbolas — r| = and the narrow bands ^ yi (inside- 
outside cascade), the support of the functions p(ri, ...,rp) in 1 -|- 1 space-time (C,'?") 
is restricted as follows: 

p-i 

p(ri, Tp) ~ 5{^p - yp) n K^i - yi)S{Ti - Tp) (2.8) 

i=l 

With this constraint and after integrating over proper time variables r^, eq.(2.7) is 
simplified as follows: 

/p 
{d^fu] exp(i qu ■ r^i^pdVi, r±i; Mt) (2.9) 
i=i 

where Pc(yj, q±i) = Pc(?i, ■-, %), PciVh r±i; Mt) = pdyi-.yp, r±i--r±p; Mt), {dVj_j} = 
d^f±i...d'^'r±p and Mt — Z^f=i Jtiti- Imposing thermal behaviour on the transverse 
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mass distribution, we may write Pc{yi,f^±i', Mt) = exp{—^)pc{yi,r±i), and since the 
Boltzmann factor does not affect the intermittency properties in momentum space 
we absorb it in pdVi, qxi) and finally obtain: 

PciVi, q±i) = / {d'^r_ii} exp(i ^ gli ■ r_Li)pc{yu r±i) (2.10) 

i=i 

where Pc{yi,r±i) represents hadronic densities in configuration space at T = T^. 

We consider the integral representation (2.10) as a model equation, based on 
rather general assumptions. It connects, at T = T^, observable density fluctuations 
in 3D momentum space {y, q±) with corresponding density fluctuations in 3D configu- 
ration space {y, f±). It is, therefore, a suitable equation for studying multidimensional 
intermittency patterns associated to self-similar structures of the hadronic system in 
real space, at T = Tc. 



3. Criticality and Fractality 

Our basic hypothesis in this treatment is a simplified picture for the hadronization 
process, in a high energy-density collision, as a static, second-order QCD phase tran- 
sition in thermal equilibrium. It generates at T = strong density fluctuations 
resulting in the development of a fractal dimension Dp in real space according to 
the standard behaviour of thermal systems undergoing such a transition and char- 
acterized by an isothermal critical exponent 6 at T = Tc {Dp = f^) 0]. In this 
picture it is natural to assume that the hypothetical second-order QCD phase transi- 
tion manifests itself as a critical behaviour of the classical 3D system (hadronic fluid) 
specified at T = Tc by the complete set of densities {pc{yi---yp, r_Li...r±p), p = 1,2, ...} 
introduced in eq.(2.10). With this interpretation, the characteristic power law of the 
chiral phase transition at T = Tc, {qq) ~ m^^^, valid in the quark phase (m is the 
quark mass), must have its counterpart in the hadronic phase, with the same critical 
exponent 6: 

{n)-p'" (^) = ^ (^ = ^^) (3.1) 

where (n) is the average multiplicity of hadrons, Qc the grand partition function at 
T = Tc, given in terms of the densities pc{yi,'^±i), and fi the chemical potential for 
the production process, normally identified with the pion mass {fi = m.,^). 

For completeness, one can easily verify that the critical behaviour (3.1) of a clas- 
sical system in D dimensions is associated to a self-similar structure in configuration 
space corresponding to a fractal dimension Dp = j^. For this purpose, consider 
the grand partition function Qdp, V) of such a system, at T = Tc, in terms of the 
densities Pc{xi, ...,Xp) {p = 1,2...) as follows: 

Qc(/i, V^) = E / {d''x.}pc{xu...,xp) (3.2) 

where z = exp{fi/Tc) and V is the volume of the system in D dimensions. A self- 
similar structure in configuration space with fractal dimension Dp requires: 

pc{Xxi, Xxp) = X^^^-^^^pcixi, ...,Xp) (3.3) 
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and in the limit /i — > eqs.(3.2) and (3.3) lead to a scaling property, Qc{fi,V) = 
F(/iV^^^/^), as follows: 

Qcif^,v) = Y.(^^^^Y^ L {^^e:}Pc(6,...,e;) (3.4) 

p V Tc J pi J\^,\<i 

In the thermodynamic limit {V — > oo) the scaling form Qc = F^fjV^^/^) combined 
with the normal boundary condition Q^. ~ exp[/(/i)l^], leads to the critical behaviour 
(3.1) with 5 = 75^5^, giving finally Dp = This simple relation of the critical 
exponent 5 with the fractal dimension Dp oi the system at T = T,, together with 
eq.(2.10) provides a basic ingredient in our approach since it makes the connection of 
the hadronic intermittency with a second-order phase transition [|^ more transparent 
and easier to understand. 

In this context and in order to establish a link between intermittency patterns in 
momentum space with criticality of the hadronic system at T = Tc one has to search, 
using eq.(2.10), for the most general class of fractals in momentum space {y,q±) 
generated by a self-similar structure in 3D configuration space {y, rj_) corresponding 
to a fractal dimension D^ = For this purpose and in order to impose self- 

similarity, we introduce new coordinates in both spaces, appropriate to the fact that 
the system, being finite, is not necessarily translational invariant. We have: 

r±i = f±i-R± (z = l,...p-l) R± = -{f±i + ...+r±p) 

q±i' = q±i-Q±i (z = i,...p-i) Q± = -{qLi + --- + q±p) 

p 

y\ = y^-Y {i = l,...p-l) Y = \y^ + ... + yp) (3.5) 

With these new variables, the problem of finding fractals in momentum space, asso- 
ciated to self-similar structures in configuration space is formulated as follows : 

p,(y:,g-l,';y,Q^) = j d^R^e^'^^'^^ j{d^ru'} 

p-i 

■ exp(z af,„q±f, 'f±„ ')pc{y'i, r±i Y, R±) (3.6) 

PciXy'i, R^) = A(^^-3)(P-iVc(2/:, '\ Y, R^) (3.7) 

PciXy-, Y, Q^) = A(^^-3)(P-i)p,(y:, gl, Y, Qx) (3.8) 

We have used eq.(2.10) together with appropriate self-similarity conditions corre- 
sponding to fractal dimensions D3 and D^ in configuration and momentum space re- 
spectively. The coefficients a^^ in eq.(3.6) are fixed by the transformation equations 
(3.5). From eqs.(3.6) to (3.8) one easily obtains stronger self-similarity constraints as 
follows: 

p^iXy,, gl, y, Q^) = A(^«+^-^)(^-i)p,(y:, u^ Y, Q^) (3.9) 
X'qu Y, Q^) = A(^-^-2)(^-i)p,(i/:, u^ Y, Q^) (3.10) 
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or in a unified form: 

~PM., y. Q±) = A^(^3+D3-4)(p-l)^i(z53-i.3-2)(p-l) ~^(^.^ g^) (3,11) 

The geometrical interpretation of eq.(3.11) is that the fractal in "iD momentum 
space, generated by a fractal F3 in configuration space, has necessarily the structure of 
a Cartesian product of two fractals, Fi and in rapidity and transverse momentum 
space respectively, F3 = Fi x F2 [§. The fractal dimensions involved in this structure 
are: Di = ^(-Ds + -D3 — 2), D2 = ^(-Ds — -D3 + 2) satisfying the Cartesian product 
relation Di + D2 = D3. Similarly, in configuration space one obtains: 

pM,^^r±^ ';Y,R^) = Ai(^^+^^-^)(^'-^)/x^(^^-^^~2)(P-^)p,(l/:,f^, ';Y,R^) (3.12) 

showing that the fractal F3 has also the structure of a Cartesian product, F3 = 
FiX F2, of two fractals, one in rapidity {Fi) and the other in impact parameter space 
(F2). The corresponding fractal dimensions in this case are: Di = ^{D^ + — 2) , 

D2 = \{D^ -D3 + 2) and D3 = Di + D2. 

The realization of this geometry in terms of particle densities in momentum space 
leads to the general form: 

pM, qu Y, Q^) = pMi)pc2{q±i ')m{y', Y, Q^) (3.13) 

~ (3) 
where pci, Pc2 obey self-similarity and is invariant, under scale transformations : 

y- -> Xy-, q±i fiq±i as follows: 

pMv:) = A(^-^)(^-^Vci(i/0 

Pc2ipq±^') = /i(^^-')(^-')pc2(gl. ') 

<l>^l{Xylpq±^■,Y,Q^) = <^^fJMq±^■,Y,Q±) (3.14) 

In summary, we have shown that the most general fractal solution of eq.(2.10) in 
configuration and momentum space, associated to a second-order phase transition 
during the development of the collision, has the structure of a Cartesian product 
in both spaces. The geometrical characteristics of these fractals are summarized as 
follows: (a) in configuration space, F3 = Fi x F2 with = Di + D2 and (b) in 
momentum space, F3 = Fi x F2 with Fi = Fi, Di = Di, D2 = 2 — D2, D3 = 
2 + Di — D2- It is of interest to note that the whole pattern of fractals in both spaces 
depends on two universal indices {Di, D2) which are constrained by the isothermal 
critical exponent 5, Di + D2 = 

In momentum space, the pattern of fractals has to be completed by considering 
ID projections of F2 within its support space (transverse momentum plane). Physi- 
cally it is equivalent to search for ID intermittency patterns in azimuthal angle (50) 
or transverse momentum {Sq±) domains. According to deterministic fractal geometry 
01 if D2 < 1 the projection of F2 onto ID subspaces of the transverse momentum 
plane is also a fractal F^i with the same fractal dimension D±i = D2. Otherwise, 
if D2 > 1, the ID projections of F2 do not develop any fractal structure {D±i = 1). 
This property remains valid for random fractals as well and can be verified for F2 as 
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follows: 

Consider the projection of F2 onto any ID w-space by integrating the densities 
Pc2{<li.i ') over the orthogonal f -space (m, v are usual Cartesian coordinates in 2D 
transverse momentum). The projected densities are given by the integral: 

P±liui) = J Pc2{uuVi)dvi...dvp-i (3.15) 

which, after imposing the self-similarity constraint (3.14) on pc2, gives: 

pW(A«.) = A(^-^)(-^)pe(«.) (3.16) 

For < 1 eq.(3.16) represents a self-similar structure in ID w-space with fractal 
dimension D±i = D2 whereas if > 1 the singularity for A ^ disappears and 
the projected particle densities have a smooth non-fractal behaviour. We therefore 
conclude that there are, in general, two possible classes of fractal patterns in momen- 
tum space, depending on the actual value of the fractal dimension D2, which can be 
attributed to a critical behaviour of the hadronic system at T = Tc. Although D2 > 1 
cannot be excluded on geometrical grounds alone, it can be easily seen that it requires 
rather unrealistic values for the critical exponent S and, therefore, can be safely ig- 
nored on physical grounds. In fact, from the general relations D2 = 2 — D3 + Di and 
Ds = (0 < -Di < 1), one obtains the bound D2 < which in the sector D2 > I 
implies 6 < 2. Since it is very unlikely for the critical exponent 6 to get such small val- 
ues in any realistic 3D system undergoing a second-order phase transition one may 
consider solution with D2 < 1 as the only acceptable class of fractals in momentum 
space associated to a second-order quark-hadron phase transition. The geometrical 
characteristics (embedding dimensions, Cartesian products, fractal dimensions Dp) 
of this solution are summarized as follows: 



3D 


F,{y)xF2{q^,(j)) 


Dp 


= 2 + D,-D2 


2D 




Dp 


= 2-D2 


2D 


F±i{q±) X h{y),F^i{<P) X h{y) 


Dp 


= 2 + ^1-^2 


ID 


Fxi(g±),F±i(0) 


Dp 


= 2-D2 


ID 




Dp 


= D, 



The physical content of the geometrical pattern (3.17) may be revealed by studying 
the behaviour of factorial moments in small domains of momentum space, using ap- 
propriate self- similarity relations implied by the above fractal structures. A set of 
power laws, dependent on two universal parameters {Di,D2), is expected to be gen- 
erated, resulting in a multidimensional intermittency pattern with a linear spectrum 
of indices. The details of this treatment are given in section 4. 



4. Multidimensional Intermittency 



The scaled factorial moments, Fp{6^l), in small domains of momentum space are 
suitable for probing geometrical structures of form (3.17) and therefore connecting a 
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second-order phase transition at the level of quark-gluon dynamics with observable 
multiplicity fluctuations at the level of the hadronic S'-matrix. The basic ingredients 
in this approach are the integrals: 

(51]3) ~ m^r^ [ ~p,{y[, qu Q^){dy[}{<fqu '}dYSQ^ (4.1) 

{p = 2,3...) where integration is restricted in the vicinity of a fixed configuration: 
yl = 0, q±i ' = 0, Y = Yq Q± = Q±q {i = 1, 2...,p — 1). In this domain eq.(4.1) is 
simplified as follows: 

(5^3) ~ ms)'-^ I PM. q±i Yo, Q±o){dy'^{d'ui '} (4-2) 

The moments F^^^S^l^) refer to 3D intermittency but similar quantities F^^^ and 
Fp^^-* in one or two dimensions are easily obtained by projecting out onto appro- 
priate momentum subspaces. Using the multiplicative form (3.13) for the densities 
PciUi, Q±i '] Q±o) we may write a multidimensional pattern of moments as follows : 

F^'\6n,) ~ {6Q,Y-^ I pMd^M yo){dy[] (4.3) 

F^'\m - mf-' I , pS(«.)$Jnl(«.; (4.4) 
F^'\5n',) ~ my-'l^, Pci(2/DpS(^.)*Snl(i/:,«.;>^o,Qf ){rfy:}{rfu.} (4.5) 

(5r]2) ~ {5n,y-^ f p,,{q^, ')^fl{qu Q±o){d'q^^ '} (4.6) 

F('\sn,) ~ {6n,y-p f pMi)~Pc2{qu ')^^M. q±^ ^o, Q±o){dy:}{d'u^ '} (4.7) 

where Q^^ is the component of Q±q along the w-direction in the transverse momentum 
plane. The integration volumes SQd are identified as follows: SQi = Sy , Sfl[ = Sq± 
or (50, 6Q2 = Sq±6y or 6(f)6y, 6Q2 = S(j)6q± and 6Q3 = 6y6(f)6q±. Their size is 
taken 6Qd ^ with M ^ 1 and the dependence of moments on M reveals the 

intermittency effects in the process under consideration. 

The integrands in eqs.(4.3) to (4.7) obey self-similarity corresponding to codi- 
mensions D — Dp = 1 — Di, D2 — 1, D2 — Di, D2,l + D2 — Di. Therefore using the 
transformations y'^ = M~^9i,q±i ' = M~^Si,Ui = M~'^'ji, we obtain the power laws 

Fi'\5n[) ~ M(^-^)(-^)/,p2(7.)$«(7.;gr){rf7.} 

Fi'\5n'2) ~ M(^-^^)(*^-i) / p,i(^^,)pS(7,,)«fS.(^„7^;>?.,Qf){^^.}{^7.} 

Fi'\sn2) ~ M^^(^-i)/ Pc2im?JM;Q^o){d's,} 

J SuJ2 

Fp(3)(5fi3) ~ M^i+^^-^^^^^'-i) 

■ / Pcli0^)Pc2im^nM. ^f, Yq, Q ±o){de^} {d^ S,} (4.8) 
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where the integration domains have finite, M-independent volumes Suoj:, ~ M^SVL^. 
In summary, we may write the complete multidimensional intermittency pattern cor- 
responding to fractal structures (3.17) as follows: 

\aF^^\5y) = (1 - Di)(p - 1) InM 

In F^^\6(f)) = lnF«(5gx) = (^2 - l)(p - 1) InM 

In Fl,'\6q^,6y) = In fI;'\6<I), 6y) = {D 2 - D,){p - 1) In M 

\nFp{5(P,5q^) = D2{p-l)lnM 

\nFf\5y,5(l),5q^) = {I + D2 - D^)ip ~ l)lnM (4.9) 

with the constraint Di + D2 = {S > 2). A number of characteristic properties 
of the pattern (4.9) can be easily observed, independently of the actual values of the 
universal indices Di,D2: (a) There is a certain degree of degeneracy at the level of ID 
and 2D moments, given by the equalities Fj^^^{6q^) = F^^\5(t)) and F^^\5q^,5y) = 
F^'^\6(f), 6y) in the limit M > 1 (5?/ ~ M-\ 6q_L ~ M-\6(f) ~ M~^), (b) introducing 
the index t] = as a measure of intermittency strength, eqs.(4.9) with b > 2 

show that the weakest effect appears in \D intermittency in rapidity (// = 1 — -Di) 
and the strongest in 2D intermittency in transverse momentum plane (r/ = and 
(c) the observed universality of ZD intermittency index at present experim [T^ may 
be attributed to the universal character of Di,D2 and in Fialkowski's notation, one 
may write 02 = |(1 + -D2 — Di). 

In summary, we have found, under a certain number of simplified assumptions, a 
complete multidimensional intermittency pattern dependent on two universal indices 
Di,D2, which, we claim, is the most general pattern associated to a second-order 
quark-hadron phase transition. 



5. Phenomenology 

At present experiments the energy density in hadron-hadron or nucleus-nucleus col- 
lisions is not high enough to generate events of hadronized quark-gluon plasma in 
abundance and therefore the universal pattern (4.9) has to be modified by adding 



non- universal terms within the context of a two-component model ||T3[ . These terms 
represent conventional, non-singu correlations which contribute to the factorial mo- 
ments as follows |[13[| : 

)(M) = a4^)m(^-^-)(^-i) + )(M) (5.1) 

where G^^\M) are slowly varying functions of M, and Ac <^ 1 is the probability for 
producing quark-gluon plasma in the actual experiment (Ac = ^^)- In a simplified 
version of this model one may assume that the conventional components G^^\M) 
are constant (M-independent) terms G^^^ and thus obtain a generalized Fialkowski's 
pattern as follows [p^j : 

ln(Fp(^) - )) = {P - 1) In(M^) + constant (5.2) 
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valid for all moments {p = 2,3,...) and all dimensions {D = 1,2,3). Hence the 
modified linear forms (5.2) are suitable for phenomenological comparison with present 
experiments in order to reveal the universal indices t] = of the pattern (4.9). 

In the particular case of 3D intermittency, Fialkowski has verified the universality of 
the index r/ by comparing the linear form (5.2) for p = 2, D = 3 with measurements 
corresponding to a number of processes. In his notation (r/ = 02 for p = 2 and 
D = 3) the universal index 02 is fixed by the data in the range 02 = 0.4 — 0.5 W^ - 
For phenomenological purposes we may now express the fractal dimensions (Di, D2) 
in terms of the universal indices (02,5), using the relations: Di + D2 = f 
i(D2-^i + l). We find: 



35 



-<5' 
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Equations (5.3) together with the constraints < Di < 1, 1 < D2 < 2, 6 > 2 put a 
restriction on the index 02: 

-3 + rM 3-rT5 ^'-'^ 

leading to a maximal range of values allowed for 02, | < 02 < 1. The lower bound 
gives a measure of 3D intermittency strength and it is remarkable that the phe- 
nomenological values of 02 lie in this domain, close to the lower bound |]12|. The 



bounds (5.4) on the other hand show a qualitative interrelation of the measurable 
index 02 with the critical exponent 5 and may lead to a phenomenological estimate of 
the allowed values for this fundamental parameter. In fact, for two distinct cases (a) 
6 = 3 corresponding to the mean-field approximation and (b) 5 ~ 5 corresponding 
to the universality class of the n = 4 Heisenberg magnet (M ~ Hs) suggested by 
Wilczek [0, we have (a) < 02 < y| and (b) | < 02 < | respectively. We observe 
that Fialkowski's phenomenological estimate of 02 is consistent with the mean-field 
approximation value 6 = 3 but marginally consistent with the n = 4 magnet value 
(5 ~ 5. It seems that present data favours rather small values of the critical exponent 
6 < A, close to the mean-field approximation value 6 = 3 [Fig.l]. For completeness we 



Figure 1: The bounds eq.(5.4) in a 5-02 plot together with Fialkowski's values for 02 
(shaded region). 

may write similar bounds for ID and 2D intermittency indices, rj = ■ For this 
purpose and according to the pattern (4.9) we put 771 = \ ~D\ {6y), rii± = D2 — I 



or 6q±) for ID intermittency and 772 = {6y6q±_ or 6y6(t)) r]2± = ^ {6q±6(l)) for 

2D intermittency. Following the same procedure as in the case of 3D intermittency 
we find, for a given value of the critical exponent 6, the bounds 

36 36 

ID : 0<,n< 3-—^ ; ^-2 < < 1 
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showing a strong 2D intermittency effect in the transverse momentum plane (| < 
'72± < !)• Finally, the phenomenological capacity of our approach is illustrated in 
Fig. 2 considering as a basic ingredient Fialkowski's universality for 3D intermittency 
at the level of the second factorial moment F2 (M). In Fig. 2a the consistency of 
the bounds (5.4) with the trend of the data is shown and the accumulation tendency 
towards the lower bound 02 = | is illustrated. We have used the outcome of a recent 



analysis of Fialkowski, based on NA35 measurements [|I4| related to a number of 
processes {OAu, SS, SAu). In Fig. 2b a particular multidimensional intermittency 
pattern is singled-out corresponding to the choice 5 = 3 (mean-field approximation) 
and 02 = 0.44 (average value from Fialkowski's analysis) in the general solution (4.9). 
This illustration suggests that a phenomenological exploitation of our model would 
require an extension of Fialkowski's analysis taking into consideration measurements 
of ID and 2D patterns as well. In such a treatment one could be able not only to 
reveal a universal multidimensional intermittency effect associated to a second-order 
quark-hadron phase transition, but also to measure the critical exponent 6 since, for 
a given value of 02, the slopes of ID and 2D moments in Fig. 2b are sensitive to the 
actual value of 6. 



Figure 2: (a) The moments (5.2) for p = 2, D = 3 corresponding to the limiting slopes 
02 = 1/3, 1 and (b) a particular multidimensional intermittency pattern eq.(4.9) with 
5 = 3, 02 = 0.44, are compared with Fialkowski's data analysis. The power laws are 
arbitrarily normalized at the same point. 



6. Conclusions 

In this work we have studied the general structure of multidimensional intermittency 
patterns in hadronic collisions generated by a second-order quark-hadron phase tran- 
sition. For this purpose we have searched for the most general class of fractals in 
momentum space associated to self-similar structures in configuration space nor- 
mally developed at T = Tc [0. The significance of self-similarity in configuration 
s pace (space-time) for intermittency phenomena has been recently emphasized by 
Bialas in connection with the Bose- Einstein interference effect [0]-A basic ingredient 
in our approach is the standard inside-outside cascade picture for the development 
of the collision, in the hadronization sector, allowing the use of rapidity as a lon- 
gitudinal coordinate both in configuration and momentum space (^ ~ y). In this 
framework, employing the reduction formalism for the inclusive ^'-matrix, an integral 
representation connecting hadronic densities in momentum and configuration space 
was established. We have found that the general fractal solution both in momentum 
and configuration space has the geometrical structure of a Cartesian product of two 
fractals, one in rapidity space and the other in the transverse-momentum or impact- 
parameter space respectively. The resulting multidimensional intermittency pattern 
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has a linear spectrum of indices, as expected, dependent on two universal parame- 
ters {Di,D2), geometrically identified with the fractal dimensions in rapidity (-Di) 
and impact parameter space {D2). Physically these universal indices are constrained 
by the isothermal critical exponent 6 of the strongly interacting system at T = Tc, 
Di + D2 = For realistic values 5 > 2 we have seen that ID projections of the 
fractal system in the transverse-momentum plane are also fractals with fractal dimen- 
sion D±i = 2 — D2- Hence ID intermittency in azimuthal angle (50) or transverse 
momentum {6q±) reflects the development of a 2D fractal structure in transverse mo- 
mentum plane with fractal dimension D2 = D±i < 1. In this case ID intermittency 
(in 6(j) or Sq±) is weaker than the corresponding effect in 2D {S(j)Sq±) in accordance 
with previous proposals about the dependence of intermittency strength on dimen- 
sionality in momentum space [0. On the contrary, ID intermittency in rapidity is 
genuine effect since the fractal structure Fi in the Cartesian product F = Fi x F2 
eq.(3.17) is an intrinsic property of the critical system and not the artefact of a struc- 
ture in higher dimensions. Its strength, in particular, depends on the basic index Di 
and is not generated by projecting out fractals embedded in higher dimensions. 

Our approach to intermittency, based on the critical properties of the strongly in- 
teracting system, is suitable for hadronic collisions (hh) and especially for relativistic 
heavy-ion processes {AA) with high energy- density. It is complementary to pertur- 
bative QCD-inspired treatments, mainly applicable to e'^e" intermittency where a 
branching mechanism at the parton level is expected to dominate the density fluctu- 



ations of the produced hadrons [|T^. 



The phenomenological implications of our solution eq.(4.9) can only be revealed 
after subtracting the conventional contributions to factorial moments in the frame- 
work of a two-component model |]13|. A simplified version of this treatment leads 



naturally to Fialkowski's universality in 3D; in order to derive a direct connection 
of measurable intermittency indices with the critical exponent 5 at T = Tc we have 
argued that an extension of Fialkowski's analysis in one and two dimensions is needed. 

Finally, it is of interest to note that a natural extension of our approach from static 
to dynamic scaling hypothesis could easily accommodate non-trivial modifications of 
the intermittency pattern eq.(4.9), at least at the level of second moments F2^\ due 



to the presence of the dynamic critical exponent |]18l . It is tempting to associate such 
modifications with a more realistic picture of inside-outside hadronization process as 
a critical phenomenon out of equilibrium 0, [l^ 



In conclusion, we have shown that multidimensional intermittency phenomena in 
hadronic collisions may originate from a simple pattern of fractal structures (mono- 
fractals) in configuration space (space-time) characteristic of a second-order phase 
transition. A careful study of these phenomena is therefore necessary, both theoreti- 
cally and experimentally, in an attempt to establish a link between criticality of the 
underlying theory (QCD at T = T^) and complexity of the hadronic system at the 
transition temperature (S"- matrix at T = Tc). 
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